Analyzing the spin-bath model without simulations 
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On the basis of a lemma designed to decide whether a discrete system decoheres or not with no 
need of computer simulations, in this paper we analyze the well-known spin-bath model. The lemma 
allows us to predict the decoherence of the system by analytical means. 



I. INTRODUCTION 



In previous papers we have developed a general theoretical framework for decoherence ([l[, Q, Q, Q), which can be 
applied to open and closed systems. The conceptually relevant step in this framework is the selection of the relevant 
observables, relative to which the question about decoherence is posed. Then, decoherence can be explained in three 
steps: 



O 
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1. First step: The space Or of relevant observables is defined. 

2. Second step: The expectation value (Or) p ^, for any Or 6 Or, is obtained. 

3. Third step: It is proved that (0 R ) p u\ — (0 R ) pR n) reaches a final equilibrium value: 
OT lim (0 R ) p{t) = lim (0 R ) Mt) = (0 R ) P , = (O r ) pr , MO R e O r (1) 

-I— > 

As it is explained in another paper of this issue (@), this general framework strictly applies when the limit of 
eq. ([T]) exists, and this happens when the Riemann-Lebesgue theorem is valid. But the Riemann-Lebesgue theorem 
LJ 1 strictly applies only in cases of continuous energy spectrum. Nevertheless, it can be shown that the validity conditions 
of the discrete analogue of the Riemann-Lebesgue theorem are expressed by the following lemma: 

Lemma 1 Let {xt} be a set of points uniformly distributed, and f{x{) : R — > K be a discrete function defined over 
■^j- ' {xi}, such that: 

in ' 

(N ■ • i e [0,N] and N > 1. 

cn, 

• 3G e N, 3P e N such that P > 1 and {xt} = \J k=1 {x^-i)^p+i)+i, —■> x k(P+i)} = Uk=i^k- In this 
will say that the set {xi} is quasi- continuous of class 1. 

• VXk, f(x rk ) = Ck, with x rk £ Xk- In this case we will say that f(xi) G C 
then, 

N 

' <=0 

where tp is the recurrence or Poincare time. 

The general aim of this paper is to apply this lemma to the spin-bath model, in order to show that this method allows 
us to predict decoherence by analytical means. For this purpose, the paper is organized as follows. In Section 2 we 
will explain how the general framework for decoherence applies to open systems, in particular, to models traditionally 
treated by means of the environment-induced decoherence (EID) approach In Section 3 we will present the 

spin-bath model, showing which the relevant observables are in this case and how the expectation values have to be 
computed. Section 4 will be devoted to compare two methods for analyzing the model: the standard method, based 
on computer simulations, and the analytical method, based on our lemma. Finally, in Section 5 we will draw our 
conclusions. 
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II. EID FROM THE GENERAL FRAMEWORK 



In the case of the EID approach, the three steps of the general framework for decoherence are usually not explicit 
in the formalism. However, the theory can be rephrased in such a way that it can be analyzed from that framework. 
In this section we will undertake this task in order to apply our just introduced lemma in a following section. 

1. First step: Let us consider a closed system U that can be decomposed into a proper system S and its 
environment E. Let the Hilbert space of U be H — Hs® He, where Hs is the Hilbert space of S and He the Hilbert 
space of E. The corresponding von Neumann-Liouville space of U is C — H ®H — Cs® Ce, where Cs = Hs® Hs 
and Ce = He® He- A generic observable belonging to C reads 

O = °s } ® °e ] G A with { P G C s and G C E (3) 

i.e. O is an observable with coordinates {0\^-g) — O^jO^l), where ... are the indices corresponding to Hs, 
and a,/3, ... are the indices corresponding to He- The relevant observables are those having the following form: 

Or = Os ® Ie G Or, with coordinates (Oij5 a p) (4) 

where Ie is the identity operator in Therefore, 0^ C C is the subspace of the relevant observables, in this EID 
case those essentially corresponding to the proper system S. 

2. Second step: The expectation value of any observable Or G Or in the state p of U reads 

{Or) p = Tr (pO R ) = Piajf) °» <W = X] °H PiajP 5 <*P = °H Pi<*j<* ( 5 ) 
ijafi ij a/3 ij a 

The reduced density operator pr is defined by tracing over the environmental degrees of freedom, 

Ps = Tr E peC' s , with coordinates ^2 Pi a ja\ = (Pij) (6) 

where C' s is the dual space of Cs- Therefore, the expectation value {Or) p u) can be expressed as 

(Or) pW = Tr (p(t) Or) = Tr (p(t)(O s ® I E )) = Tr (p s (t) O s ) = (O s ) Ps{t) (7) 



3. Third step: The EID approach studies the time evolution of the reduced density operator ps(t) governed by 
an effective master equation. For many physical models where the space Or has a finite number of dimensions, this 
approach shows that, for t —¥ oo, ps(t) reaches an equilibrium state ps*'- 

Ps(t) —> ps* (8) 

Since ps* is obviously diagonal in its eigenbasis, the system S decoheres in the eigenbasis of ps*, which turns out to 
be the final decoherence basis. But if we take into account the definition of ps as a partial trace (see eq. (O), we can 
obtain the limit of the expectation values of eq. ([7} as 

lim {O s ) Ps{t) - Urn (0 R ) p(t) = (O s ) PSt = {0 R ) P , (9) 

where p* is such that ps* results from the projection of p* onto Or. Therefore, for any observable Or C Or, 

hm (0 R ) m = (Or) p , (10) 

This result can also be expressed as a weak limit 

W - lim p(t) = p* (11) 

t— >OG 

If the just obtained eq. (fTU|) is compared with eq. (TTJ) , it turns out to be clear that the EID approach can also be 
formulated from the viewpoint of the closed composite system U and, from this perspective, it can be explained in 
the context of the general framework introduced in the Introduction. In other words, the split of the closed system 
into a proper open system and an environment is just a way of selecting the relevant observables of the closed system. 

The limit of eqs. © and (JSJ) can be computed by means of two different strategies: 
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. By solving the unitary evolution equation for p(t), computing Tr (p(t) Op) and finding the limit. 

. By solving the non-unitary evolution equation for ps(t) = Tte (p(t)), computing Tr {ps{t) Os) and finding the limit. 

Of course, the two strategies give the same result: although the second is the usual method in the EID literature, 
the first may lead to a simpler solution, as in the model developed in the next section. 



III. THE SPIN-BATH MODEL FROM THE GENERAL FRAMEWORK 



The spin-bath model is a very simple model that has been exactly solved in previous papers (see [8] ) . Here we will 
study it from the general framework applied to the EID approach, as presented in the previous section. This task will 
allow us to compare the method traditionally used in the literature for solving the model with the method based on 
our lemma. 

Let us consider a closed system U = PUP\ U. . . UPjv = PU (U^LjPj), where (i) P is a spin-1/2 particle represented 
in the Hilbert space Hp, and (ii) each Pi is a spin-1/2 particle represented in its Hilbert space Hi. The Hilbert space 
of the composite system U is, then, 



N 



H = H P ' 



(12) 



In the particle P, the two eigenstates of the spin operator S P ^f in direction v are l^) , |JJ»): 

1 



In each particle Pi, the two eigenstates of the corresponding spin operator in direction v are |ti) , \U) 



Si,!? |fi) = ~ It,} 

Therefore, a pure initial state of U reads 

\i>o} = (a\ir) + b\itt) 



\U) = \U) 



N 



)(ai\ti)+/3i\U)) 



(13) 



(14) 



(15) 



where I a 1 2 + |6| 2 



1 and \a.i 



1 (for a generalization with M spins {ff, JJ,} and N spins {t,|}, see g). If 



the self-Hamiltonians Hp of P and Hi of Pj are taken to be zero, and there is no interaction among the Pi, then the 
total Hamiltonian H of the composite system U is given by the interaction between the particle P and each particle 
P (see 1, 0): 



1 N 



gi(\ti)(ti\-\U)(U\) 




(16) 



where Ij = + is the identity operator on the subspace Hj. Under the action of H, the state jf/'o) 

evolves into 



|^(t))=o|^>|^(t)>+6|^>l^(t)> (17) 

where 

N 

Mt)) = M-t)) - <S> ("« e ~ mt/2 ItO + ^ e i9it/2 II,)) (18) 
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A. Computing the expectation values 

The space O of the observables of the composite system U can be obtained as O = Op <g> (®f = iOi) : where Op is 
the space of the observables of the particle P and Oi is the space of the observables of the particle Pj. Then, an 
observable 0^0 = %®% can be expressed as 

JY 

= P ®((g)O l ) (19) 

where 

op = s n it) (n\ + s u m + 8 U m (t\ + 8 U m <w e o P (20) 

Oi = eg 1 ti> (t< I + eg I I + 4? I W (ti I + 41 1 ti>Ui I G ^ (21) 

Since the operators Op and are Hermitian, the diagonal components s^, s^jj., eg, eg are real numbers, and the 

off-diagonal components are complex numbers satisfying = s|^, eg = eg*. Then, the expectation value of the 
observable O in the state \ip(t)) of eq. (fTT)) can be computed as 

<C% (t) = (|a| 2 Sl y t + |6| 2 SM )r (t) + 2Rc[ a fe* SM r 1 (t)] (22) 



where (see 



ro(*) - II [l^! 2 4t + Iftl'eg +2Re(a !; /3*ege^*)] (23) 
i=l 

JV 

ri(t) = II [l^| 2 4t eWt + lAI'ege"^* +2Re(a i ^eg)] (24) 



j=i 

B. Selecting the relevant observables 

In the typical situation studied by the EID approach, the open system S is the particle P, and the remaining 
particles Pi play the role of the environment E: S — P and E = ufL 1 Pi. Therefore, the relevant observables Op of 
the closed system U are those corresponding to the particle P, and they are obtained from eqs. (fl9)) , ([20]) and (|21l) . 

by making eg = eg = 1 and eg = 0: 

Op = O s ®Ie=\ s »As){s'\ J ® (25) 

\s,s'=ft,4 / \i=l / 

The expectation value of these observables in the state \^}{t)) of eq. dl~7l) is given by 

(0 R ) m = \a\ 2 s n + \b\ 2 s u + 2 Re[a&* s m r(t)] (26) 

where 

N 

r{t) = (smwrn = n (i^i 2 e ~ mt + i&i 2 e * 9!t ) ( 2? ) 

i=l 

This means that, in eq. (J22), r (£) = 1 and r x (t) = r(*). 

IV. SIMULATION VERSUS PREDICTION 

A. Simulation: the usual method 

In order to know the time-behavior of the expectation value of eq. (|21)1) . the time-behavior of r(t) has to be 
computed. From eq. (121)1) |r(i)| 2 results 

N 

Ht)\ 2 = Y[{\a l \ i +\(3 l \ i + 2\a l \ 2 \f3 l \ 2 cos2g l t) (28) 
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If \ai\ 4 and \f3i\ z are taken as random numbers in the closed interval [0, 1], such that \cti\ 2 + \/3i\ 2 = 1, then 
rnax(K| 4 + |/3»| 4 + 2\a,\ 2 \fr\ 2 cos 2g t t) = ((| a ,| 2 + |&| 2 ) 2 ) = 1 
min(|a l | 4 + |/3 I | 4 + 2|a l | 2 |/3 l | 2 cos(2 5l t)) = ((H 2 - |/3,| 2 ) 2 ) - (2 |«,| 2 - l) 2 (29) 



Therefore, (|«i| 4 + |/?i| 4 + 2|ai| 2 |/3i| 2 cos 2gtf) is a random number which, if t 5^ 0, fluctuates between 1 and 



(2|cd 2 -l >2 



In order to obtain the limit of r(t) for £ — > 00, different numerical simulations are performed and presented in 
the literature, where the aleatory numbers \(Xi\ and |/?,-|~ are obtained from a generator of aleatory numbers: the 
generator fixed the value of \oti\ , and the is computed as |/3j| = I — \oti\ . The value of the gi and the time 
interval [0,io] for the computations is usually stipulated. In general, the model is studied and the conclusions about 
decoherence are drawn by means of this kind of numerical simulations (Q, [||). 

B. Prediction: using the lemma 

The first step for applying the lemma is to express the expectation values in the energy eigenbasis. In this model the 
energy eigenbasis coincides with the eigenbasis of the spin in direction z, which is obtained in terms of the tensorial 
product of the cigenstates of the spin in direction z for all the particles. In other words, the eigenvectors of H, which 
form a basis of %, are 

It) |tl> - It*) - ItiV-l) Itjv) 

It) |ti)...|tfe)-ltiv-i)M 

|^> \U) - \U) ... \In-i) \In) (30) 



Then, it is easy to sec that in eq. (j!6[) H is written in its diagonal form, and that the expectation value of eq. 
is expressed in the energy eigenbasis. 

In order to simplify the expressions, we will introduce a particular arrangement into the set of those eigenvectors 
by calling them \Ai): the set {|.4i}} is an eigenbasis of H with 2 N+1 elements. The \Ai) will be ordered in terms of 
their eigenvalues, which depend on the number of particles of E having spin \\) in any state. Then, we have: 

• Two states where all the particles of E have spin |f): 

\A,) = |1bt,-,t,t) =>H\Ai) = \ (iZito) I A) 
|^_x) = |t, t, ti t) => H \A-i) = -\ (E£i3i) l-^-i) (31) 

• 2~N states where only one particle of E has spin ||): 

\Aj) = % t, t, I, t, t, t) H \A 3 ) = \ (YZi9i ~ 9k) \Aj) 

\A- 3 ) = \^A,.-.AAA,-.-AA)^H\A- 3 ) = -±(z?=i9*-9k) \A- 3 ) 

with j = 2,3,..., N+ 1 and k = 1,2,..., N (32) 

• (N — 1) N states where two particles have spin |T): 

\Aj) = |t, t, , t. i, t, , t, I, t, t, t) H \Aj) = l - (T,tx9i - 9k - 9i) \Aj) 

\A-j) = |t, t, , t. U t, , t, I, t, t, t) =► H \A-j) = - 1 - (E,= i9i - 9k - 91) \A- 3 ) 

(N — 1) N 

with j = N + 2,N + 3,...,N + l + ± —!— and k,l = 1,2,..., N (33) 
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• For the remaining particles with more spins \\), the procedure is analogous. 

Let us consider the two extreme cases. If the coupling coefficients gi are random numbers, in principle all eigenvectors 
are different. If the coupling coefficients are gi = g, then there are 

• , ■ , N 

2 eigenvectors with eigenvalue —g 

N - 2 

2N eigenvectors with eigenvalue — - — g 



2 — — — — eigenvectors with eigenvalue — - — g (34) 

I J. V L J • t . zj 

with I = 0, 1, ...N. In this case, H is degenerate: it has 2 N+1 eigenvectors but only 27V different eigenvalues. In both 
cases, random gi or equal gi, the number of different possible energies is large when TV is large enough. 
On the other hand, eq. (|2l)|) can be written as 

(0 R ) m = ]T PiOi + PtOu e iw "' (35) 

i v 

where 

^PiOi = \a\ 2 s n + \b\ 2 s u (36) 



N 

Y / PtO v e^ t = 2Re[ab*s M l[(\a l \ 2 e-^ t + 1(3^6^)} (37) 

v i—1 

The r.h.s. of this last expression includes a binomial product which can be rewritten as a sum by means of the 
following strategy. First we define the index v that establishes the number of the term of the sum: since in eq. (|3"7| 
there are 2 N terms, then v = 0, 1, • • • 2 N — 1. Then we define the number p v ^ as the i digit of the number v written 
in the binary system. Moreover, each term of the sum is a product of N exponentials of the form e~* 9i *, which can 
be grouped into a single exponential e luJut . The w„ are all possible additions and subtractions between the coefficients 
gi] so, a generic w„ can be computed as 

^ = (eL (-ir-ft) (38) 

Precisely 

v = 
v = 1 
v = 2 
v = i 
v = A 



v = 2 N - I = I ■ ■ ■ l b — ► uj 2N = (-EtiSi) (39) 
On the basis of this strategy, we can define the discrete function /d(u;„) as 

N 

= n i>^i 2 ( 4 °) 
fc=i 



• • • 000b 
• ■ 001b 
0---010 b 
— Ollf, 
• • • I00f, 



w 

,W3 



i=i 9i ~ 9n 



Ei=i 3i - 9jv-i + J 

V-iJV-2 \ 

E l= i & - 9n-x - 0jv J 

w 4 = ( Ei=l 3i - 3JV-2 + ffJV-l + .9w 
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where 7^ = (ak — Pk)Pv,k + fik, which is equal to au if p v ,k — 1 an d is equal to (3k if ?Vfc = 0. Then, the binomial 
product of eq. (|3"T)) results 



N 2"-l 



[] (|a 4 | 2 e- wt + |A| 2 e wt ) = ]T fd^e**"* (41) 

then, 

2 N -1 

i i/=0 

In order to apply our lemma to eq. (|42|) . it is necessary that /d(w„) £ L\. On the one hand, since {w„} has 2 N 
elements, then for N » 1 the set {w^} is quasi-continuous of class 1. On the other hand, since fd(w v ) is defined in 
eq. P0|) . where < |7i/,/c| 2 < 1, then fd{u v ) is the product of N numbers lower than 1. Therefore, if TV ^ 1, 

N 

i7.,fei 2 <i^/dK)=ni>.fei 2<<iifiv>>i o</ d K)< £ «i 
fc=i 

^ max(|/ d (w M ) - < 1 (43) 

But this is precisely the condition for fd{w u ) £ C\. As a consequence, according to the lemma, 

Y^PlO.e^^Q =>■ (O fi ). 0W ^E^ ( 44 ) 

v i 

and we can conclude, with no need of computer simulations, that the system decoheres. 



V. CONCLUSIONS 



In another paper of this issue, a discrete analogue of the Riemann-Lebesgue theorem is presented and, on this basis, 
a lemma relevant for discrete models is introduced: such a lemma provides a criterion for deciding whether or not 
the system decoheres with no need of numerical simulations. In order to present an example of how these results can 
be usefully exploited for the study of decoherence in discrete models, in this paper we have applied that lemma to 
the well-known spin-bath model, and we have shown that the conclusion drawn from that application agrees with the 
results obtained by means of computer simulations in the previous literature. 
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